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Abstract. Conditions within the Homunculus nebula around Eta Car are de¬ 
termined by many of the same physical processes that occur in molecular clouds 
in the interstellar medium. But there is one major exception - we know when 
the ejection occurred and something about its composition and initial state. 
The gas was warm, ionized, and dust-free when it was located within the star’s 
atmosphere and it is currently cold, molecular, and dusty. It undertook this 
transformation in a bit over 150 years. It offers a laboratory for the study of 
physical processes in a well-constrained environment. We derive a photoion¬ 
ization model of the Homunculus nebula that reproduces many of its observed 
properties. We conclude by outlining how observations of the Homunculus could 
address basic problems in the physics of the interstellar medium. 


1. Introduction 

Eta Carinae, one of the most luminous stars in the Galaxy, offers a laboratory in 
which a variety of physical phenomena can be studied. The star has undergone 
at least one episode of substantial mass loss, and is likely to end its life as a 
supernova. Understanding the physics that occurs within the ejecta will offer 
insight into the initial stages of chemical enrichment of the interstellar medium. 
Both molecules and grains are known to have formed in the Homunculus nebula, 
gas that was ejected in the nineteenth century. The physical conditions in this 
nebula are the subject of this paper, along with a sketch of how Eta Car might 
be used as a test bed to understand critical interstellar processes. 


2. An illustrative model of the Homunculus 

The Homunculus nebula is seen in the optical by reflected starlight and is one of 
the brightest objects in the sky in the mid-IR, showing that substantial amounts 
of dust are present (Westphal & Neugebauer 1969). In the infrared, emission 
lines of [Fe ii] and H 2 are seen and have been traced across the nebula in long- 
slit observations (Smith 2002). These show a double-shell structure, with Fe"*" 
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present in the inner parts of the Homunculus, and H 2 in the outer zone. Addi¬ 
tionally, the Homunculus shows a double-shell structure in the dust color tem¬ 
perature distribution, with a cool outer shell at ~140 K and a warmer inner 
shell at ~200-250 K (Smith et al. 2003). These observations help constrain any 
model of the envelope. 

We model the shell as a constant-density layer with an inner radius of 
1.7x10^^ cm and a thickness of 10^^ cm, appropriate for material in the wall of 
the southeast polar lobe along our line of sight to the star (Smith 2002; Smith et 
al. 2003). The electron density in the inner Fe"*" region, deduced from infrared 
[Fe ii] line ratios, is ~10^ cm“^ (Smith 2002). The hydrogen density must be 
substantially higher here since the gas is mainly neutral. From the geometry of 
the nebula and the presumed total mass of ~12 Mq for the entire Homunculus 
(Smith et al. 2003), we infer a hydrogen density of 10^ cm Of course, this 
depends on the gas-to-dust mass ratio, taken to be 100 (but see below). For 
these parameters the column density through the shell is A^//=10^^ cm“^. 

The dust that is clearly present in the nebula will be the catalyst in forming 
the observed H 2 . The measured visual extinction is Ay ~ 4 mag, (e.g., Davidson 
& Humphreys 1997), corresponding to an extinction per unit column density of 
Ay/Nu ~ 4x 10“^^ mag cm^. The extinction is observed to be grey and the dust 
color temperature in the shell is near the equilibrium blackbody temperature 
(Whitelock et al. 1983; Hillier et al. 2000; Smith et al. 1998, 2003), suggesting 
that the grains are large. We assume that the grains are similar to those seen in 
Orion’s Veil, R=5, and assume that only the silicate component is present. This 
large R is produced by a grain size distribution that is lacking in small grains, 
which will affect the H 2 formation rate and grain photoelectric heating. 

The ratio of C/0 is less than unity in the ejecta, (Davidson et al. 1986), 
suggesting that the chemistry will be dominated by oxygen-bearing species once 
formation of CO is complete, and that the chemistry will eventually lead to 
oxygen-rich solids, motivating our use of the silicate grain type. The observations 
of silicate features in the infrared (e.g., Gehrz et al. 1973), and the absence 
of a graphite feature in the ultraviolet (Viotti et al. 1989) supports this idea. 
For simplicity we leave the silicate dust to gas ratio at its ISM value, which 
corresponds to an extended source Ay/Nn ~ 9 x 10“^^ mag cm^. The grain 
size distribution and dust to gas ratio will affect the details of our calculations, 
as well as clumping of the material, but not the overall results. 

The assumed gas-phase abundances are listed in Table 1, along with other 
parameters. Relative to H, He is overabundant, while O is highly underabun- 
dant, presumably due to partial CNO cycling. The N/H ratio, roughly ten times 
solar, corresponds to the conversion of nearly all C and O into N (Davidson et 
al. 1986; Smith & Morse 2004). Most C is expected to be in the form of the ^^C 
isotope. 

We assume that the stellar continuum is represented by an interpolated 
20,000 K CoStar atmosphere with a total luminosity of SxIO^Lq. We add a 
high-energy component corresponding to a 3x 10® K blackbody with a luminosity 
of 30 Lq (e.g., Corcoran et al. 2001). The lack of a prominent H ii region shows 
that few hydrogen-ionizing photons strike the inner edge of the nebula, most 
likely due to photoelectric absorption by the stellar wind. We extinguish the 
net continuum by photoelectric absorption due to a neutral layer of 10^^ cm“^ 
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to account for this. Some high-energy photons are transmitted and they help 
drive the chemistry. The incident stellar continuum is shown in Figure 1. 

We also include the galactic background cosmic ray ionization rate. The 
actual ionization rate may be higher if radiative nuclei are present. Cosmic rays 
have effects that are similar to X-rays - they provide ionization that helps drive 
the chemistry. 


3. Calculations 

We simulate the conditions within the nebula using the development version of 
Cloudy, last described by Ferland et ah (1998). 

Recent updates to Cloudy include an improved molecular network that 
allow for calculations deep in molecular clouds. Some of these improvements are 
discussed in Abel et al. (2004). The developmental version of Cloudy currently 
predicts molecular abundances for ~70 molecules involving H, He, C, O, N, Si, 
and S. Approximately 1000 reactions are in the network, with most reaction 
rates taken from the latest version of the UMIST database. Our predictions are 
in good agreement with other codes that are designed to predict conditions in 
PDRs. 

The computed ionization, thermal, and molecular structure are shown in 
Figures 2 and 3. The emitted continuum is the solid line in Figure 1. 

Hydrogen is predominantly atomic at the illuminated face of the cloud. We 
assume that no H-ionizing radiation escapes from the stellar wind, but the H"*" 
density at the illuminated face is quite sensitive to the transfer of the incident 
stellar continuum in the Lyman lines. If it is bright in these lines then hydrogen 
can become ionized by a two-step process. An excited state is populated by 
absorption of a Lyman line, which then decays into the H® 2s level. The Balmer 
continuum can photoionize atoms in this state, creating a thin region with H'*' 
and a high electron density. The Lyman lines quickly become self-shielding and 
the process is no longer important, although a small amount of H+ is produced 
across the nebula by cosmic ray and X-ray ionization. 

H 2 forms at depth of ~3.5xl0^® cm, where the Lyman-Werner bands be¬ 
come optically thick, the continuum between La and the Lyman limit is heavily 
extinguished, and the destruction rate of H 2 goes down dramatically. As Fig¬ 
ure 1 shows, little light escapes at short wavelengths. Grains are the dominant 
opacity across the cloud, helping shield H 2 and allowing for efficient formation 
by catalysis on grain surfaces. The Fe^ profile is also shown, indicating an anti¬ 
correlation between Fe"*" and H 2 . Observations (Smith 2002) show that Fe^ and 
H 2 are segregated, occupying the inner and outer zones of the Homunculus walls, 
respectively, which generally agrees with the structure in Figure 2. 

The formation of large amounts of H 2 initiates the formation of heavy- 
element molecules (see Figure 2 and Table 2). H 2 is a step in the formation of 
H^ and Hg", the highly reactive ion-molecules that undergo ion-neutral reactions 
to form molecules containing heavier elements. Large amounts of CO form when 
its electronic bands become self-shielding. For this calculation, CO fully forms 
at depth of ~9xl0^® cm. We assume a C/0 abundance ratio of 0.5. At the 
shielded face the n{CO)/n{Ctot) ratio is nearly unity and a significant amount 
of O is in the form of OH. 
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Nitrogen is strongly enhanced in the ejecta, and several nitrogen-bearing 
molecules are shown in Table 2 and Figure 2. As expected, N 2 is the domi¬ 
nant molecule, although only a small amount of N is in this form - N remains 
predominantly atomic. 

The gas temperature is shown in Figure 3. It lies in the range ~50 K to 
~300 K and is typical of a PDR. The temperature is mainly maintained by 
a balance between grain photoionization heating and cooling by fine structure 
lines of C, O, Si, and Fe. The temperature falls at the point where H 2 forms 
due to the strong absorption by its electronic transitions and also by atomic C. 
In the coldest regions heating by cosmic rays becomes important, together with 
cooling by CO rotation transitions. 

The Homunculus appears “lumpy”, showing that the envelope does not 
have constant density or pressure. These calculations show two possible sources 
of local instabilities which might help form blobs. The radiative acceleration, 
mainly caused by the absorption of the incident continuum by grains and lines, 
exceeds 3x10“^ cm s“^ at the illuminated face but falls to 3x10“® cm s“^ at 
the shielded face. This suggests that significant radiative acceleration can occur 
over the ~10^ yr lifetime of the ejecta. It may be Rayleigh-Taylor unstable 
because of the decreasing acceleration. The thermal balance is a second source 
of instability - the temperature derivative of the net cooling, defined as cooling 
minus heating, is negative across much of the ejecta; this material is thermally 
unstable. Both instabilities will be the focus of future work. 

The emitted spectrum is shown as the solid line in Figure 1. Thermal 
infrared emission from grains is prominent, along with CO lines in the mm. The 
10 ^m silicate feature is present since the calculation only included a silicate 
component. A detailed comparison with observations would help constrain the 
model further, and so help deduce properties such as the dust abundance and 
composition. 


4. The Future 

Our main purpose is to point to directions for future work on the Homunculus, 
with the goal of using it as a laboratory to understand basic physical processes. 
The Homunculus is an especially well-posed problem. The ejecta were once part 
of a hot stellar atmosphere and so must have initially been warm, ionized, and 
dust-free. Today it is cold, dusty, and molecular. How did it go between these 
states during the time since its ejection? 

This dust could not have formed in the atmosphere of Eta Car itself, since 
the energy density temperature is above the condensation temperature of most 
solids. Thus the dust now seen in reflection is most likely to have formed in 
the material after being expelled from the star. The cycle of dust destruction 
and formation is still poorly understood (Draine 1990). These newly formed 
dust grains have a large size, and measurements of the extinction curve across 
the spectrum would help quantify their radii. A comparison between intensities 
of the thermal IR continuum and H 2 or CO lines would measure the dust-to¬ 
gas ratio. Infrared spectral features can also reveal the dust composition. This 
is especially important in light of recent work suggesting that supernovae are 
important sources of new grains in the galaxy (Morgan et al. 2003). 
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The Homunculus is predicted to be predominantly molecular. A molecular 
inventory could be obtained from UV observations of electronic absorption lines 
or from the many prominent rotational emission lines that are expected in the 
IR - mm. We have a good idea of the initial gas-phase chemical composition 
so this inventory would test current chemical reaction networks, and especially 
the theory of H 2 formation on grain surfaces. The chemistry will be affected by 
additional cosmic ray ionization if radioactive nuclei are present in the ejecta, 
and also by the isotopic variations - most C should be rather than Can 
the chemistry test these assumptions? 

Finally, the prominent structures seen in the Homunculus can test dynam¬ 
ical theories. Could thermal or radiation driving instabilities play a role? 
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TABLE 1 - Model parameters 


Parameter 

(Homunculus) 


He/H 


C/H 


N/H 


0/H 


Remaining 

elements 


riH 


Inner radius 


Thickness 


Filling factor 


Covering factor 


Grains 


Cosmic ray 



1 


0.03 


R=5 Orion silicate 


Galactic 

Background 


Parameter 

(Stellar) 


Wind extinction 


Atmosphere 


Jkf _ 


L 


Corona 



Atmosphere 


Jkf _ 


L 


TABLE 2 - Predicted column densities (cm 


Species 


Log N 

Species 

Log N 

Species 

Log N 

22.53 

N 2 

14.64 


12.55 

22.52 

CH 2 

13.98 

HCO+ 

12.50 

20.18 

CH 

13.45 

NO 

12.48 

18.86 

CH3+ 

13.27 

SiO 

12.42 

18.50 

HCS+ 

13.09 

HCN 

12.41 

18.27 

cs 

13.06 

HeH+ 

12.23 

18.24 

CN 

12.88 

SiN 

12.21 

18.17 

CH3 

12.87 

OH 

12.18 

18.13 

CH 4 

12.80 

^S^ 

12.12 

17.54 

H 20 

12.63 

"CHj 

12.02 
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Figure 1. The incident stellar continuum, emitted by the central star, is 
shown as a dashed line. The predicted emergent continuum is the solid line. 
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Figure 2. The predicted molecular, atomic, and ionic densities vs depth into 
the shell. The upper panel shows that the inner half of the shell emits Fe II 
while H 2 is found in the outer two-thirds, in agreement with observations. 

The lower panel gives densities of the major molecular constituents. In the 
outer-most region CO is the dominant molecule, while N 2 is the dominant H- 
carrying species. The text identifies the remaining molecules. At the shielded 
face the molecules, in order of decreasing abundance, are CO, N 2 , CH 2 , H 2 O, 

SiO, CS, NO, H+, CH^, SiN , HS+, and CH+. 
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Figure 3. Figure 3, The predicted electron temperature and density struc¬ 
ture. The total hydrogen density was assumed to be constant. 
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